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Abstract 

To compare two space-times on large domains, and in particular 
the global structure of their manifolds, requires using identical frames 
of reference and associated coordinate conditions. In this paper we use 
and compare two classes of time-like congruences and corresponding 
adapted coordinates: the harmonic and quo-harmonic classes. Besides 
the intrinsic definition and some of their intrinsic properties and dif- 
ferences we consider with some detail their differences at the level of 
the linearized approximation of the field equations. The hard part of 
this paper is an explicit and general determination of the harmonic 
and quo-harmonic coordinates adapted to the stationary character of 
three well-know metrics, Schwarzschild's, Curzon's and Kerr's, to or- 
der five of their asymptotic expansions. It also contains some relevant 
remarks on such problems as defining the multipoles of vacuum solu- 
tions or matching interior and exterior solutions. 
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1 Introduction 



Let us consider three well-known vacuum solutions of Einstein's equations: 

1. The Schwarzschild solution written in Droste-Hilbert coordinates: 

dsl = -fi-—) de + \ df' + f\de' + sin2 ed^'), (1) 

with 

f > 2M. (2) 

2. The Curzon solution written in spherical Weyl-related coordinates: 

dsl= - e-'^^'/'df 



exp I — ) [dr + r do ) + r sm 6 dip 



(3) 



with 

f > 0. (4) 
3. The Kerr solution written in Boyer-Lindquist coordinates, 
, 2 / 2Mf\ , 2 4Mr 



= — 1 — dt H —a sin 6' (it d(p 

+ ^rf^2^^2^^2^ ^2^^2^ 2Mr^2gi^^2]gi^^2^^2^ 

with 

= f2 + a2 cos > 2Mf, A = + - 2Mf > 0. (6) 

We have used different notations to distinguish the radial and angular 
coordinates to emphasize that there is no connection a priori that makes 
sense across the three solutions. 

These three space-time metrics share some intrinsic properties, i.e. which 
are independent of the system of coordinates being used to describe them. 
Namely: the three solutions possess a main global time-like Killing vector 
in the domains specified in Eqs. (@), (§) and (||). This vector is actually a 
generator of a group of isometrics which includes as sub-group the group of 
rotations around an axis, and the specular symmetry across a plane. The 
intrinsic differences are also well known. Namely: the Schwarzschild and 
Curzon solutions are static, i.e. ^" is integrable, while the Kerr solution is 
only stationary. The Schwarzschild solution is spherically symmetric while 
the Curzon and the Kerr solutions are only axially symmetric. 
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The three systems of coordinates used in Eqs. (Q), and (||) share also 
some common properties: 

i) The space coordinates are adapted to the main time-hke Kilhng vector 
^" as well as to the space-like Killing vector corresponding to the axial 
symmetry. This is reflected by the fact that the gravitational potentials are 
independent of t and ip. They are also adapted to the specular symmetry of 
the solutions. 

ii) They clearly show their Euclidean behavior at spatial inflnity in the 
frame of reference corresponding to the congruence deflned by 

These common properties by no means make these coordinates unique, 
and therefore it would be even heuristically unjustifled to refer to these three 
metrics as being written using a "common" system of coordinates. This paper 
deals with the problem of further restricting the systems of coordinates up 
to the point of making such an assertion acceptable. 

We shall see that the solution to this problem sheds new light on the 
problem of matching them to interior solutions. On the other hand it is 
obvious that being able to use a "common" system of coordinates is the 
only available means to compare different solutions on large domains of their 
manifolds. 

Coordinate speciflcations can be made on different criteria. The Droste- 
Hilbert coordinates in Eq. (|I|) are completely determined by the requirement 

dsl\^=f\de^ + sm^ed^^), (7) 

S being any of the 2-surfaces t = const., f = const. It can be used only for 
spherically symmetric space-times. The Weyl-like coordinates in Eq. (|]) are 
specially tailored for static axially symmetric metrics. The Schwarzschild 
solution could be written using such coordinates but we all know how mis- 
leading this would be from a physical point of view. 

The Boyer-Lindquist coordinates are akin to Droste-Hilbert coordinates 
in the sense that when the parameter a in the Kerr solution is made zero 
then Eq. (j^) becomes Eq. ([l|). 

It is clear that the connections between the three systems of coordinates 
we have considered dbTQ clS yet too loose and partial to be of any interest. 
Above all because they do not shed any new light on any other problem. 

We take in this paper the point of view according to which to determine 
the physical content of the space-time metric being considered as well as to be 
able to compare two of them on large domains it is necessary to use systems 
of coordinates whose deflnition make sense independently of the metric to 
which the deflnition is applied. 

In Section we consider the deflnition and some properties of harmonic 
and quo-harmonic congruences and the intrinsically related harmonic and 
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quo-harmonic coordinates. The hnear approximation discussed in Section ^ 
provides a simphfied framework to illustrate some interesting properties of 
those congruences. 

Section ^ deals with the problem of writing the Schwarzschild solution, 
up to terms of order included, both in harmonic and quo-harmonic co- 
ordinates adapted to the intrinsic symmetry structure described above. 

In Section ^ we consider the same problem for the Curzon solution. Our 
results show in particular that Schwarzschild and Curzon solutions coincide 
when terms of order and higher can be neglected. This is the case when 
M is the mass of the Sun and the physical system being considered is the 
solar system. 

The same problem for the Kerr solution is analysed in Section ^. Our 
results show in particular that when considering asymptotic developments in 
1/r of this metric using quo-harmonic coordinates it is necessary to include 
some logarithmic terms behaving as r~^lnr. 

2 Harmonic and quo-harmonic coordinates 

Whenever the metric 

ds^ = gaf3 dx'^dx'^ (8) 

of a space-time model is given using an explicit system of coordinates x", two 
associated geometric structures are also implicitly given: i) the congruence 
C of world-lines defined by the parametric equations 

X* = constant, i = 1,2,3, (9) 

which we shall always assume to be time-like in the domain V of interest, 
and ii) the foliation JF defined by the equation 

= constant, (10) 

which we shall always assume to be space-like on T>. 

Two questions can be asked which are relevant in the theory of frames of 
reference or more generally whenever we are interested in the suitability of the 
coordinate system to handle a particular problem: i) does the congruence C 
belong to any particular type, like being for instance a Killing or Born or any 
other type intrinsically defined? and ii) can the foliation be intrinsically 
characterized independently of C, or does it have any particular intrinsic 
properties connected with it? 

Systems of coordinates can be divided in two classes: those that are spe- 
cific to particular models or classes of models and are useful because they take 
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advantage of particular properties of them, and those that make sense what- 
ever model is being considered. Isotropic, Droste-Hilbert or Schwarzschild 
coordinates for space-times with spherical symmetry, or Weyl coordinates for 
static axially symmetric models belong to the first class. 

Gauss coordinates belong to the second class, but to our knowledge the 
single problem where they have played an important role in the development 
of the theory of General Relativity is in the problem of matching exterior 
to interior solutions of Einstein's equations across space-like hypersurfaces. 



More on that on Section 4.3 



2.1 Harmonic congruences and coordinates 

A widely used system of coordinates belonging to this second class is 
that of harmonic coordinates, which is characterized by these two groups of 
equations: 

□x* = (11) 

and 

□x° = 0, (12) 

where □ is the intrinsic d'Alembertian of the space-time metric. The splitting 
above in two groups of the four equations Dx" = will be understood in a 
moment. 

Conversely, given a time-like congruence defined by a unit tangent vector 
field it may be asked what type of adapted coordinates can be most 
appropriate to use. For instance, we know that if is coUinear with a 
Killing vector field then it is always possible to use a system of harmonic 
coordinates such that the congruence with parametric equations (|^) is the 
same as that defined by -u". But in general this is not the case, because if 
we consider the harmonic coordinates x* as functions of generic coordinates, 
= riy")) then the functions p must satisfy the system of equations 

□ f = 0, M"5,f = 0, rank(9,f) = 3, (13) 

and this requires integrability conditions on the vector field that in gen- 
eral are not satisfied. A particularly interesting case in which they do not 
hold is that of irrotational Born (or rigid) congruences when they are not 
Killing |jl|. As a matter of fact, Eqs. (|l^) can be considered as the definition 
of a new intrinsic type of congruences: that of harmonic congruences, i.e. 
those admitting a set of adapted harmonic coordinates of space. 

Notice that the time coordinate can always be required to satisfy Eq. ([T2|) 
and therefore this condition is, without any supplementary conditions, unre- 
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lated to any particular congruence: it is enough to find a solution of 

af = 0, ^ 0. (14) 

It is for this reason that we preferred to refer to harmonic coordinates split- 
ting the conditions on the space and time coordinates. 

2.2 Quo- harmonic congruences and coordinates 

As we have already mentioned these considerations are particularly relevant 
in any theory of frames of reference that requires to extend this concept be- 
yond the Born congruences which are, generically, notoriously exceptional. 
As a contribution towards an appropriate generalization it has been intro- 
duced H, ^ a new type of congruences which are defined much in the 
same way as we defined harmonic congruences by Eqs. ([TB|), except for a 
slight modification of the first group of conditions, which becomes 

□r+A°a«r = o, (15) 

where 

A" = -u^Vpu" (16) 

is, up to the sign, the intrinsic curvature of the congruence u". We shall also 
consider later a quo-harmonic time coordinate defined by 

af + A-dJ' = 0, ^ 0. (17) 

This class of quo-harmonic congruences contains, in contradistinction to 
the class of harmonic congruences, the whole class of Born congruences. And 
it provides sufficient generality to be an essential ingredient to define the con- 
cept of rigidity without which no useful meaning can be given to the concept 
of frame of reference. This is why we shall spend some time to justify the 
consideration of quo-harmonic congruences and quo-harmonic coordinates; 
two concepts that are related but that by no means are identical. 

2.3 Harmonic and quo-harmonic coordinate classes 

The next section is dedicated to compare Einstein's linearized equations using 
harmonic or quo-harmonic coordinates. In the remaining sections instead we 
shall restrict ourselves to the consideration of the three stationary metrics 
mentioned in the Introduction. These metrics are stationary and have been 
written in a system of coordinates adapted to a main Killing vector. All 
Killing congruences are both harmonic and quo-harmonic. This can be seen 
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as follows. In a system of coordinates adapted to a Killing vector field all 
the potential ga/3 are independent of time and Eqs. (|15D become 



A/^ = -^9, i^^jgg''dj-') =0, g^ det , (18) 

where it has been used the decomposition 

ds"^ = - ^ dt + ipi dx"^ ^ + gij dx'dx^ (19) 



with 

i = V-goo, V^i = C^9oi, 9ij = 9ij + ^Vi'/'i- (20) 

Since Eqs. (|18D always have |^ three solutions independent of t it follows 
that the Killing congruences are quo-harmonic. On the other hand Eqs. (p!3| ) 
become 



Af = —d, (v^ rd,f) = 0, gi, ^ e9^,, 9 = det (g,,) , (21) 

and these equation also have three independent solutions independent of t. 
Therefore Killing congruences are also harmonic. 

Notice however that if Aj 7^ harmonic coordinates are different from 
the quo-harmonic ones even for a congruence that is both harmonic and quo- 
harmonic that includes as we have seen the Killing congruences — in 
which case they reveal different aspects of the three models being considered. 

The harmonic coordinates are derived from those used in Eqs. (|I]), (^) 
and (H) by a pure space transformation 

x' = x'{y^,X,fi), (22) 

where is any of the systems of coordinates used in ([I|), (^) and (^, A is 
the set of parameters on which depends the metric (M, but also a for the 
Kerr metric) and /i is any set of constants coming from integrating the co- 
ordinate conditions being demanded. For to be harmonic and Cartesian 
at space infinity the gravitational potentials ga/sit, x*) have to be solutions of 
the equations 

r^^ ^ ^7"^r^^ = 0, (23) 

which are Eqs. ([TT|) when harmonic coordinates are used. Here F^^ are the 
Christoffel symbols. Of course any other system of coordinates derived 
from x* by 

= f\x^) (24) 
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with functions /* independent of the parameters A and n can be said to belong 
to the harmonic class and be used to fulfill the requirement of universality 
mentioned before. In particular, we shall systematically use not harmonic 
Cartesian-like coordinates but rather polar coordinates that belong to the 
same harmonic class and are related to them by the familiar formulas: 

= r sin 9 cos ip, = r sin 9 sin ip, x^ = r cos 9. (25) 

We will also use quo-harmonic coordinates, which are derived from the 
original ones by space-like transformations of the type ( p2D but requiring 
instead of (^3]) Eqs. (pTsD, which become 

= + g'^^Aj = (26) 

when quo-harmonic coordinates are used. For stationary metrics these equa- 
tions reduce to 

T'' - g^^dj In ^ = 0. (27) 

3 The linear approximation 

We assume in this section that the metric admits a congruence and a 
system of adapted coordinates such that it can be written as 

gal3 = Val3 + ha/S, (28) 

where rjap are the Galilean coefficients of the Minkowski metric, and h^js are 
small quantities whose products can be neglected, as well as the products of 
derivatives of any order, in the domain T> of interest. 

Adapted coordinate transformations leaving form-invariant (^) are up 
to a Lorentz transformation 

= x° + C°(a;"), x" = x' + C{x^), (29) 

where the are small quantities, and the C do not depend on time. 

Under such an adapted coordinate transformation the quantities ha/s 
transform as follows: 

h'oo = /ioo + 2aoCo, (30) 
h'oi = hoi + diCo, (31) 
h'ij = hij + diCj + djCi. (32) 
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The differential invariants under these transformations are 



Ai = -dih^Q - dohoi, (33) 

Qij = dihoj - djhoi, (34) 

^ij = dohij, (35) 

Rijki = —-{dikhji + djihik — diihjk — djkhii). (36) 



These expressions give in this approximation the values of intrinsically well 
defined objects associated to the time-like congruence P, ^: the sign-reversed 
acceleration (|16|) , the vorticity field, the deformation rate and the Zel'manov- 
Cattaneo tensor ||. 

3.1 Harmonic congruences 

Let us consider the quantities defined in (|23|) , which in this approximation 
reduce to ^ 

Ti = -dohoi + djhl - -di{-hoo + h), (37) 

where 

h = (38) 

Under an adapted coordinate transformation those quantities and their de- 
rivatives with respect to time transform as follows: 

r,, = Ti + AO, doTi, = doTi. (39) 

Since at this approximation a congruence is harmonic iff a system of coordi- 
nates exists such that Fj' = 0, it follows that 

doT^ = - + doAj = 0, S = (40) 

characterizes the harmonic congruences to this approximation. This is in 
fact an invariant condition that guarantees the existence of a solution of the 
equation 

AC = -r. (41) 

for (i not depending of 

If the harmonic congruence is a Born congruence then Hij = and it 
follows from (^OD that Aj must be independent of time. If moreover we assume 
that Qij is also independent of time then it is easy to show that the congruence 
is in fact a Killing congruence, i.e. a system of adapted coordinates exist such 
that doga(3 = 0. 
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Let us now consider the linearized Einstein's field equations, 

Rap = Ual3, Ual3 = T^/B " ^{v'"'Tpa)'nal3, (42) 

where Raf3 is the linearized Ricci tensor: 

R^f^ = -^(D V - ^/^r^ - d^Tf,). (43) 

Unlike in Section here □ is the d'Alembertian associated to the Minkowski 
metric, T^^ is some approximation to the energy-momentum tensor, Fj are 
the quantities defined in (0), and we get from Eq. (pSf): 

To = -dohoo + djhi - ^9o(-/ioo + h). (44) 

Assuming that the congruence is harmonic, that harmonic coordinates 
are used, Fj = 0, and that a foliation with Fq = is selected, then the familiar 
equations 

Rap = -TpKp = Uap (45) 

are obtained. These show explicitly the hyperbolic type of Einstein's lin- 
earized field equations when space-time harmonic coordinates can be used. 

3.2 Quo-harmonic congruences 

Let us consider now the quantities defined in (p^): 

f, = 0,14 - ^dih. (46) 

Under an adapted coordinate transformation these quantities and their de- 
rivatives with respect to time transform as 

Ti, = Ti + AO, doTe = doTi. (47) 

Since at this approximation a congruence is quo-harmonic iff a system of 
coordinates exist such that Fj/ = 0, the condition 

dot, = d.T.) - U,T. = (48) 

characterizes the quo-harmonic congruences to this approximation. This is 
an invariant condition that guarantees the existence of a solution of the 
equation 

AO = -f, (49) 
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with Q independent of x^. 

From (^O]) and (0) it follows that a congruence can be harmonic and quo- 
harmonic at the same time iffAi is independent of time. Notice however that 
this does not mean that in this case harmonic and quo-harmonic coordinates 
are the same. 

Let us assume now that the congruence is quo-harmonic, that quo- 
harmonic coordinates are used, Fj = 0, and that a foliation is selected — as 
it can always be done by using (pT]) — such that Fq = 0, where 



f = d,hl (50) 

The later hypothesis is equivalent to assume that x° also satisfies the quo- 
harmonicity condition (|1^. Under these assumptions, Einstein's field equa- 
tions become 

Roo = ~ (A/ioo + dooh}j = Uqo, (51) 

Ro^ = -\ (^Ahoi + ^doih}j = f/o„ (52) 

Rij = {ahij + do{dihoj + djhoi) - dijhoo) = Uij, (53) 

where A is the Laplacian constructed with the 3-dimensional Euclidean 6ij 
metric. Introducing the traceless tensor 

kij = hij -^hSij, (54) 
the last group of equations splits in two groups: the scalar equation 



R= --{Oh- Ahoo) = U, R = 6'mij, U = 6'^Uij, (55) 
2 



and the tensor equation 

Rij -^RSij = Uij -^U6ij. (56) 
Taking into account Eqs. ([51|)-([53|), the scalar Eq. (^) reduces to 

- ^Ah = Uoo + U (57) 
and the tensor Eq. (BBI) becomes 



1 11- 

- {nkij + doidihoj + djhoi) - dijhoo) - -A/ioo^ij = Uij - -U6ij. (58) 
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Assuming Tap^xP) known everywhere and for all times, and satisfying the 
conservation equations 

do^T^ = 0, (59) 

h can now be obtained integrating the elliptic equation (^). After that /iqo 
can be obtained integrating the elliptic equation (0). And then h^i can be 
obtained integrating the elliptic equation (^21). Finally the single part of 
the metric that necessitates to integrate hyperbolic equations is the traceless 
piece kij appearing in Eqs (^81). We can say then that neither /iqo, nor /loi, 



nor h are propagating quantities in quo-harmonic coordinates. Only the 
wave-part of the metric may propagate. When solving the aforementioned 
equations one have to make sure that conditions Fq = and Fj = are 
satisfied, but they are compatible conditions because of conservation law (|59|) 
and the fact that from the field equations (|5T|)-(|53D one gets 

- ^Afo = d^T^ (60) 

and 

-]pt^ = ^o^Tt + \^o^to. (61) 

3.3 Stationary solutions 

If the frame of reference in which Eqs. ( ^5l ) or Eqs. (^)-(|53|) have been 
written is that corresponding to a Killing congruence implying that h^ii are 
time independent then the preceding equations can be split in two groups. 
The first group is common to both harmonic and quo-harmonic coordinates: 

-^A/ioo = ^00, (62) 

-^A/io. = U^i, 9,/ij, = 0. (63) 

The second group is 

- ^A/i,, = Uij, djivi - h)i (-/loo + /i) = 0, (64) 
when using harmonic coordinates and 

- ^ (A/ijj - dijh^Q) = Uij, djhl - ^dih = 0, (65) 



hij = hij of (|65| ) are related as follows 



when using quo-harmonic coordinates. The solutions hij = hij of ( |64D and 

hij = hij + diCj + djCi, (66) 
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Q being a solution of the equation 



ACi = ^d.hoo. (67) 

If Uai3 is known, smooth and decreasing fast enough at infinity, or compact 
with discontinuities of the first kind across the border of the support, then 
the solution to the preceding field equations are: 



and 

As it is well known the solutions thus obtained are at least of class C^. 
This is one of the reasons, among many others, including the analogy with 
electromagnetism, to require in general that the potential ga/3 be of class 
across the surface of discontinuities of the first kind in the theory of relativity. 
Surprisingly this requirement is not respected by many authors because of 
what appears to be a confusion between two connected but different concepts: 
that of metrics that can be matched and that of metrics that have been 
matched. 

The multipole structure of the potentials /iqo and /loi can be determined 
as it is standard in electromagnetic theory. The multipole structures of hij 
and hij are in general different but related by Eq. (0). Notice also that 
the multipole structure of a particular solution is independent of any system 
of adapted and admissible coordinates, harmonic, quo-harmonic or else, if 



instead of the potentials we consider the invariant quantities (p3|)-(|36|). Fur- 
thermore, as a consequence of the field equations ([6^)-(|65D for stationary 
metrics in vacuum, (^) is equivalent to 

Rij = —-dijhoQ (70) 

both in harmonic and in quo-harmonic coordinates, so that the invariant 
multipole structure in these two coordinate classes for this kind of metrics is 
equal and determined by the common values of /iqo and hoi- In the last section 
of this paper we shall come back to the linear approximation in a different 
context connected with the problem of defining the multipole structure of a 
vacuum solution of Einstein's equations when the source is not known. 
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4 Schwarzschild metric 



To write the Schwarzschild solution we are going to use as basis the following 
forms 

= dt, 00^ = dR, uo"^ = RdQ, oo^ = R sin 9 # (71) 
when we use harmonic coordinates and 



° - dt, uj^ = dr, u? = r dO, uj^ = r sin 6 dip (72) 



with quo-harmonic coordinates, so that the metric coefficients Qfj^y are defined 
as follows: 

ds"" = g^uLU^Lu". (73) 
These functions reduce to the Minkowskian values when Rot r go to infinity, 

lim g^,^ = r]f,^, (74) 

ft or r^oo 

and they will be written as 

gtMu = rii^u + h^^,, (75) 
where the h^j^i, are functions of R or r. 

4.1 Schwarzschild metric in harmonic coordinates 

We can obtain a system of harmonic coordinates preserving the spherical 
symmetry if we perform a coordinate transformation in the form (f, ^, ip) 
{R, 6, 0), where R = R{f), Q = 9, (j) = ip. In fact, the Cartesian coordinates 
associated to the set of spherical coordinates {R, G, 0) are harmonic and 
limfi^f^ R = r is satisfied if one chooses Pl 



R{f) = f-M + C 



^1 A 2M , ,^ 
In 1 — + M 



(76) 



2 V f 

where C is an integration constant, which will be discussed below in Sec- 



tion 4.3 



We can invert the coordinate transformation in series of M/R and the 
functions /i^i, of the non-null metric components written in harmonic coor- 
dinates to order five are 

M M2 Q + 2CM^ 6 + 

- ^r''^W^'^r^~^~r^^^^W 
m q-acm^ q-iocm^ 10 - 28c 

hn - 2^ + 2^ + ^—:^+ 2 'W^ 5 'W' 
M 2 2 2 
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4.2 Schwarzschild metric in quo-harmonic coordinates 

We can also obtain a system of quo-harmonic coordinates preserving the 
spherical symmetry if we perform a coordinate transformation in the form 
(f, 9, 0) (r, 6, ip), with r = r(f), 6 = 6,ip = ip. The Cartesian coordinates 
associated to (r, 6', ip) are quo-harmonic and the condition limf^oo r = f is 
satisfied if [|10|, 0] 



r(r) 



1 

2f 



(7J 



where A is an integration constant to be discussed later. Notice the remark- 
able fact that for A = the interval 2M < f < cxd corresponds to < r < cxd. 
No such thing is possible in harmonic coordinates. 

If we invert the coordinate transformation in series of M/r, the func- 
tions h^i, of the non-null metric components in the quo-harmonic coordinates 
are to order five 

- 2^ 29-2AM^ 



(-11 



M M2 1 - 2A M3 a 5 - 6A 
2 I I . 



M 9 M2 1 - A M3 
h22 = hs = 3— + -— + — ^. (79) 

4.3 Matching with an interior solution 

To analyse the meaning of the constants C and A appearing in ([761) -([77|) 
and (|78D -([79|) respectively, let us consider the problem of constructing a full 
model for the gravitational field of a spherically symmetric star. We can use 
the Schwarzschild metric (|l]) for the exterior field but we also need an interior 
solution and then both metrics have to be matched at the star surface. It 
is well known that the continuity across the matching surface of the first 
and second fundamental forms is enough to guarantee that the matching 
can be done without a surface layer [|T2|. But from our point of view the 
problem is not completely solved with this: to actually make the matching 
we also need to found for both metrics a set of "admissible coordinates" 



in the sense of Lichnerowicz in which all the derivatives of the metric 
coefficients are continuous across the star surface. Furthermore, we want 
this common set of coordinates to have the same physical meaning for the 
interior and exterior solutions: they must be members of a class of systems 
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of coordinates defined with independence of the particular problem under 
study. In particular we are interested here in the cases in which the common 
coordinates are harmonic (or quo-harmonic) for both the interior and the 
exterior metric. 

In fact, Quan-Hui Liu considered a particular interior with uniform 
density and found that the matching with the Schwarzschild metric (|l|) in 
harmonic admissible coordinates is only possible if the constant C in ([76|) 
has a given value depending only on the radius of the star surface. 



A similar result arises in quo-harmonic coordinates [Tl|]: an interior so- 



lution with constant density can be matched in quo-harmonic admissible 
coordinates only if the constant A in ([78|) has a particular value determined 
by the interior solution, the matching radius and the constant M. 

These two particular examples suggest that the integration constants that 
appear when solving the differential conditions of definition for harmonic or 
quo-harmonic coordinates will be fixed (hopefully in a unique way) only when 
the full problem of finding a common set of such coordinates for both the 
internal and external gravitational fields is addressed. 



5 Curzon metric 

We are going to write the Curzon metric in the same bases ([ri])-(|7^) we have 
used for the Schwarzschild metric. The metric coefficients g^p in Eq. ( [75| ) 
satisfy the asymptotic property ( [74D and we will use again notation (|75D, but 
now the h^^ quantities are functions of R and B or of r and 9. Moreover, 
they have the following structure: 

V = E hfJP2U = Y: hffPl (80) 

where h^^^ are only functions of or r and P„ and are the Legendre 
polynomials and associated functions in cos 6 or cos 6. 

5.1 Curzon metric in harmonic coordinates 

One may obtain a system of harmonic coordinates preserving the axial sym- 
metry and the reflection symmetry with respect to the plane z = by using 
an appropriate coordinate transformation in the form {r,6,ip) (i?, B,0), 
where R = R{f, 6),Q = Q{f,6),(j) = (p. By solving the harmonicity equation 
that the Cartesian coordinates associated to the set of spherical coordinates 
{R, 9, 0) must satisfy, and demanding that lim^^oo R = r,we obtain to order 
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five in the expansion parameter M/R: 



R 
f 



1 + 



sin^^M^ 2Ci + 2C2 - 3C2 sin^ ^ M3 4-5sin2^ . ^~M' 



6 



+ 



24 



sin^^- 



^4 



240 



24OC4 + (8C1 + 2OOC2 - 42OC3 - 84OC4) sin^ 9 



-(24OC2 - 525C3 - 6OOC4) sin^ 9 



(81) 



cos© 
cos^ 



l_lsin2^~^-::isin^^ 
2 f2 2 



C2_,__2n^^ 4-llsin2^ . ^~M^ 



24 



sin^^- 



^4 



1 

80 



I6C1 + 8OC2 - I4OC3 - 8OC4 
-(I2OC2 - I75C3 - 2OOC4) sin^ 9 



sin' 9 



(82) 



where Ci, C2, C3,C4 are integration constants. 

In this system of harmonic coordinates the functions h^^j) for the non-null 
components of the Curzon metric are, to order five, 



(0) 
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M M2 6 + 2CiM^ 6 + 4CiM5 
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(83) 



11 



M M2 6-4CiM3 e-lOCiM^ 10 - 28Ci 

2 u 2 I - I - I - 

R R^ 3 R^ 3 R^ 5 R^' 
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2 + 4C2 4 + IOC2 
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(84) 
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^ 9 'W 
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140 - 168Ci - 2IOC2 + 42OC3 - 36OC4 
^ 315 'W' 

- 56 ^^^^ 

,(0) _ KP 2C1 + C2M3 I-6C1-6C2M4 



'''22 



R R"^ 3 9 

40 - 96Ci - I2OC2 + IO5C3 



180 i?5 ' 

2 + 2C2 5 + I8C2 



3 9 
160 - 1176Ci - 336OC2 + 1365C3 - 4320C4 
1260 



,(4) 38 + 3I5C3 + 36OC4M5 

- 105 ^^^^ 



M M2 2C1-C2M3 l + 6Ci-6C2M^ 



R R^ 3 i?3 ' 9 ^4 

40 + 48Ci - I2OC2 + IO5C3 



180 R^ ' 

(2) _ _ 2 _ 7 + 6C2 _ 80 + I68C1 + 336C2 - 735C3 
~ ^3^ 9 'W 252 

- 105:R^- ^^^^ 

5.2 Curzon metric in quo-harmonic coordinates 

To obtain a system of quo-harmonic coordinates preserving the axial and 
reflection symmetries we will perform a coordinate transformation in the 
form (f, 9, (f) — >• (r, 9, (p), with r = r(f, 9),9 = 9{r, 9),ip = (p. By demanding 
limj:_>oo r = f and that the Cartesian coordinates associated to the spherical 
(r, 9, (p) satisfy the quo-harmonicity equation we obtain to order five in M/r 

r _ ^ I M sin^ 9 1 + ^1 - 2^2 - (3 - 3^2) sin^ 9 

f 2 f 2 f'^ 4 T"^ 

_ 3 + 3Ai - 6A2 - (13 - 9A2) sin^ 9 + 5 sin^ 9 

24 
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240 



9 + 9Ai - I8A2 + 6A4 

+ (59 - 6^1 + 327^2 - 288^3 - 21^4) sin^ 9 

-(115 + 360^2 - 360^3 - I5A4) sin^ — , 



(88) 



cos 9 
cos 9 



= 1 



sin^ 9M'^ 2 - 3/1, 



-3r + 



sin' 9— 



7-llsin2^ 



24 



240 



64 - 36yli + 387^2 - 288^3 - 6^4 



-(20 + 540^2 - 360^3 - 15^4) sin' ^ sin'^— , (89) 



sin'^ 



where Ai, A2, A3 and A4 are integration constants. 

In this system of quo-harmonic coordinates the functions /i^^j-* for the 
non-null components of Curzon metric are written as follows: 
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'^00 



M M2 9M3 29- 2^1 99 - 18^ 

2 3 I - I i . 
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105 r5 ' 
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(90) 
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12 280 r5 ' ^ ^ 
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4 



M 9 M2 1 - 3Ai + 3^2 19 - 27^2 M' 



4 



28 - 9Ai + 18^2 - 36^3 



90 r 



5 



(2) _ I + 3A2M3 13- 27^2 

94 - 441^1 + 882/I2 - 468^3 + 54^4 
630 

(0) _ 5-3A1-3A2M3 19- 27^2 

"-33 — 3 h - — H - H — 

r 4 r"^ D r-^ 3d 

28 - 9Ai + I8A2 - 36^3 
90 

^(2) _ 7 - 27^2 214 - 63Ai + 126^2 - 252^3 

4? ^ . (94) 

35 

Since now we have two sets of coordinates that can be considered common 
to the Schwarzschild and Curzon metrics, we can compare the two metrics. 
From the results in Section ^ we see that both metrics coincide if terms of 
order / B? (or M^/r^) and higher are neglected. This means that in most 
practical cases (such as in relativistic celestial mechanics) these two met- 
rics are indistinguishable if used in harmonic or quo-harmonic coordinates. 
Notice that if one naively identifies the coordinates used in (P and (|^) the 
difference between the Schwarzschild and Curzon metrics starts with terms 
proportional to M/f = M/r. These remarks have no special utility except 
that they demystify the apparent weirdness of the Curzon solution. 



6 Kerr metric 

To write the Kerr metric we will use again (|7TD-(|73D and the metric coeffi- 
cients g^i, will satisfy (^). The h^^ quantities of ([7^) are functions of R and 
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or of r and 6 and may be written as 

V = EC^^2^ /M2 = EM?A'^ /^03 = Etf3^'^^2m, (95) 
I I I 

where h^^'^ are only functions of i? or r and P„ and are the Legendre 
polynomials and associated functions in cos© or cos^. 
We will use the dimensionless quantity 



a = 



M 



(96) 



to make easier the comparison with the results for the Schwarzschild metric, 
which is recovered in the limit a = 0. 

6.1 Kerr metric in harmonic coordinates 

To preserve the axial symmetry and the reflection symmetry with respect to 

the plane 2; = we perform a coordinate transformation (f, 9, (f>) —> {R, O, 0), 
with R = R{f,9),Q = Q{f,6),(f) = (f>. If lim^^^oo P = f and the Cartesian 
coordinates associated to the spherical (P, 0, 0) are harmonic, the transfor- 
mation is given to order five in (now we have two dimensionless expansion 
parameters: M/R and a/R) by 
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1 sm — r- 

f 2 f-^ 6 f-^ 

I6D1 + + (3 - 6D2) sin^ ^1 + sin^ 9 



120 



24 
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(97) 



cose q;^ y;-M 

^ = 1-— sm^6l — 

cos 6* 2 



— a 



2l + ^2 . 2n 

— — sm 9 



xM3 



2 12P)2 - 3^2 sin^ ^ 2;^M^ 

-a sm 9 — T- 

8 f4 



+ - 



a 



52 - 32D2 + 2D3 + 80^4 - 4a2p'2 



(52 + I4OD4 - a;^(5 + I2D2)) sin^ 9] sin^ 9—, (98) 
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where Di, D2, D^, D4 are integration constants. 

In this system of harmonic coordinates the functions h^^l for the non-null 
r>nmr>onents of the Kerr metric are written as follows to order five: 

M M2 M3 6 + 2DiM4 6 + 4Di-2q;2m5 
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6.2 Kerr metric in quo-harmonic coordinates 

We can also obtain a system of quo-harmonic coordinates preserving the 
axial and specular symmetries by using a change of coordinate in the form 

{f,6,^) {r,6,Lp), with r = r{f,6),6 = 9{f,9),Lp = (p. The Cartesian 
coordinates associated to (r, 9, ip) will be quo-harmonic and lim^i^oo r — f 
will be satisfied if, up to order 1/r^, we choose 
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where we have again four integration constants: Si, S2, Bj, and Bi^. 

Notice that, unhke in all the preceding examples, here we do not have 
only powers of 1/r but also terms of the form In r/r^. 

Functions for the non-null components of Kerr metric appear as fol- 
lows in this system of quo-harmonic coordinates: 
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We may now compare the Schwarzschild, Curzon and Kerr metrics in a 
common (harmonic or quo-harmonic) set of coordinates defined with indepen- 
dence of the models being analysed. If we neglect again terms of order 1/ 
(or and higher, we see that the only difference between the two static 

cases and the Kerr metric is the term /1Q3 = 2Ma/B? (or /iq3 = 2Ma/r^). 
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7 Final comments 



The approximation scheme in Sections. ^, |^ and ^ is based on expanding an 
exact solution of the complete Einstein's equations about the infinity point 
of an appropriate radial coordinate. In the so-called linear approximation 
discussed in Section |] the exact solutions of the approximated (linearized) 
field equations were used. Nevertheless, by carefully selecting some terms 
from the expansions in Sections. ^ ^ and |^ one may recover exact solutions 
of the linearized theory. 



For instance, from the expressions in Section [4.1| for the Schwarzschild 
metric in harmonic coordinates, one may easily find the following solution of 
the linear approximation: 

2M , , 

hoo = (113) 

ho^ = 0. (114) 



The same happens (after replacing R with r) in the expressions in Section ^ 
for the Schwarzschild metric in quo-harmonic coordinates. Notice that we 
are here writing only the metric components that give the common invariant 



multipole structure, as discussed in Section |3T3. 

In the case of the expressions in Section |5]2 for the Curzon metric in 
harmonic coordinates, one may easily find the following solution of the linear 
approximation with more complex angular structure: 

2M 2M3 ^ 38 , , , 

^00 = ^-3:ffr^2(cose) + — — P4(cose), (ii5) 

hoi = 0. (116) 

Substituting r for R and 6 for B in this result we get the solution of the lin- 
earized field equations contained in the approximated solution of Section ^]2] 
for quo-harmonic coordinates. 

In the case of the Kerr metric, the solutions of the linearized theory we 
can extract from Section ^ are 

hoo - ^-2^P,(cose) + 2^P,(cose), (117) 
hoi = ho2 = 0, (118) 
hos = 4^AMcose)-l^P3Mcose) (119) 

in harmonic coordinates, as well as in quo-harmonic coordinates if one uses 
{r,e) instead of (P, 6). 
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